1 NGC-144/000169-199 
FOUR SAMPLE CALCULATION OF SENSOR ARRAY INDUCED PHASE ANGLE 



STATEMENT OF GOVERNMENT RTCHTS 
The Government of the United States of America has rights in this invention pursuant 
to Contract No. N00024-02C-6305 awarded by the U.S. Department of the Navy. 

TECHNICAL FTEi n 
The invention relates generally to signal processing and more particularly to 
demodulation of signals from fiber optic sensor arrays. 

BACKGROUND 

Fiber optic sensor arrays of a time division multiplexed ("TDM") system are often 
used to measure a change in a parameter, for example, acoustic vibration, fluid pressure 
variations, acceleration, and magnetic field intensity. The fiber optic sensor array employs a 
phase generated carrier with a period T to measure the change in the parameter at a given 
sampling rate. The fiber optic sensor array converts a phase angle associated with the 
parameter to an amplitude variation on an output pulse of light. 

The phase angle is measured through various demodulation techniques of the output 
pulse. Typical demodulation techniques employ a quadrature component Q and an in-phase 
component I of the output pulse. The quadrature component Q corresponds to a sine of the 
phase angle, and the in-phase component I corresponds to a cosine of the phase angle. An 
arctangent of the ratio Q/I is equal to the phase angle. The magnitude of the change in the 
parameter can then be calculated from the change in the phase angle. 

Calculation of the quadrature component Q and the in-phase component I requires 
multiple samples of the output pulse at specific intervals of the phase generated carrier. The 
phase generated carrier comprises a period that is significantly longer than a period of the 
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output pulse. The longer period of the phase generated carrier requires the samples to span 
several output pulses to obtain each required interval of the phase generated carrier. It is 
desirable to increase the rate of calculation of the phase angle to promote an increase in 
throughput. For a constant sampling rate, the number of required samples must be reduced to 
increase the rate of calculation of the phase angle. 

Thus, a need exists for reduced numbers of required samples for demodulation 
techniques of fiber optic sensor arrays that employ phase generated carriers. 

SUMMARY 

The invention in one embodiment encompasses a method. A sensor array employs a 
parameter to induce a time-varying phase angle (p on an optical signal that comprises a phase 
generated carrier. The phase angle <p is calculated through employment of only four samples, 
where all the four samples are based on the optical signal. 

Another embodiment of the invention encompasses an apparatus. A sensor array 
employs a parameter to induce a time-varying phase angle (p on an optical signal that 
comprises a phase generated carrier. The apparatus comprises a processor component that 
employs only four samples to calculate the phase angle (p, where all the four samples are 
based on the optical signal. 

A furflier embodiment of the invention encompasses an article. A sensor array 
employs a parameter to induce a time-varying phase angle <p on an optical signal that 
comprises a phase generated carrier. The article includes one or more computer-readable 
signal-bearing media. The article includes means in the one or more media for calculating 
the phase angle <p through employment of only four samples, where all the four samples are 
based on the optical signal. 
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DESCRIPTION OF THE DRAWINGS 

Features of exemplary implementations of the invention will become apparent from 
the description, the claims, and the accompanying drawings in which: 

FIG. 1 is a representation of one exemplary implementation of an apparatus that 
comprises one or more lasers, one or more optical switches, one or more phase modulators, 
one or more sensor arrays, one or more optical receivers, and one or more processor 
components for calculating a phase angle cp through employment of only four samples from 
an optical signal. 

FIG. 2 is a representation of a contour plot of a ratio R(M, p) for a previous design. 

FIG. 3 is a representation of one exemplary contour plot of a ratio R(M, P) for the 
implementation of FIG. 1. 

FIG. 4 is a representation of another exemplary contour plot of the ratio R(M, P) for 
the implementation of FIG. 1. 

FIG. 5 is a representation of a contour plot of the ratio R(M, P) and a total harmonic 
distortion for the previous design of FIG. 2. 

FIG. 6 is another representation of the exemplary contour plot of the ratio R(M, P) of 
FIG. 3 and a representation of a total harmonic distortion for the implementation of FIG. 3. 

FIG. 7 is another representation of the exemplary contour plot of the ratio R(M, P) of 
FIG. 4 and a representation of a total harmonic distortion for the implementation of FIG. 4. 

DETAILED DESCRIPTION 
Turning to FIG. 1, an apparatus 100 in one example comprises a plurality of 
components such as computer software and/or hardware components. A number of such 
components can be combined or divided in the apparatus 100. An exemplary component of 
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the apparatus 100 employs and/or comprises a set and/or series of computer instructions 
written in or implemented with any of a number of programming languages, as will be 
appreciated by those skilled in the art. 

Referring to FIG. 1, the apparatus 100 in one example comprises one or more lasers 
5 102, one or more optical switches 104, one or more phase modulators 106, one or more 
sensor arrays 108, one or more optical receivers 110, and one or more processor components 
1 12. In one example, the apparatus 100 demodulates an optical signal to measure a change in 
a parameter, as described herein. The laser 102 in one example comprises a continuous wave 
laser. The laser 102 generates and sends an optical signal through the optical switch 104 and 

10 the phase modulator 106 to the sensor array 108. 

The optical switch 104 in one example comprises a time division multiplexed 
("TDM") switch. The optical switch 104 gates the optical signal such that the optical signal 
comprises a stream of optical pulses. The phase modulator 106 impresses a phase generated 
carrier ("PGC") 114 on the stream of optical pulses. For example, the laser 102, the optical 

15 switch 104, and the phase modulator 106 cooperate to create one or more optical pulses 1 16 
that comprise the phase generated carrier 114, as will be understood by those skilled in the 
art. The optical pulse 116 comprises a period Tpuke- The period Tpuise in one example is 
approximately between 50 nanoseconds and 500 nanoseconds. The phase generated carrier 
1 14 in one example comprises a period Tpgc and a modulation depth of M. The period Tpgc 

20 comprises a relationship with a frequency fpgc = 1 / Tpgc, as will be understood by those 
skilled in the art. The frequency fpgc in one example is approximately between 20 KHz and 
100 KHz. The phase generated carrier 1 14 is associated with a demodulation phase offset p. 
The phase generated carrier 114 creates a time-varying phase angle equal to 
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The sensor array 108 in one example comprises one or more sensors 124, 126, and 
128, for example, mismatched path interferometers. The sensor array 108 splits the optical 
pulse 116 into one or more optical pulses 118, 120, and 122, for example, one pulse per 
sensor. The optical pulses 1 16, 1 18, 120, and 122 in one example are substantially the same. 
The sensors 124, 126, and 128 of the sensor array 108 receive the optical pulses 118, 120, 
and 122, respectively. The sensors 124, 126, and 128 of the sensor array 108 in one example 
employ one or more parameters and the optical pulses 118, 120, and 122 to create one or 
more respective interference pulses 130, 132, and 134. Exemplary parameters comprise 
acoustic vibration, fluid pressure variations, acceleration, and magnetic field intensity. For 
example, the sensor 124 splits the optical pulse 118 into a first portion and a second portion. 
The sensor 124 employs the parameter to induce a time-varying phase angle 9 on the first 
portion of the optical pulse 118, relative to the second portion of the optical pulse 118. The 
sensor 124 recombines the first portion of the optical pulse 118 with the second portion of the 
optical pulse 124 to create the interference pulse 130. A time-varying amplitude variation of 
the interference pulse 130 represents the time-varying phase angle cp between the first portion 
and the second portion of the optical pulse 118. 

The optical pulses 116 comprise an intermediary spacing such that the interference 
pulses 130, 132, and 134 comprise a relatively small spacing, for example, a high duty cycle, 
as described herein. The interference pulses 130, 132, and 134 comprise a period 
substantially equal to the period Tpuise of the optical pulse 116. The sensor array 108 sends 
the interference pulses 130, 132, and 134 to the optical receiver 1 10 in a pulse train 136, for 
example, in a serial fashion. For example, the optical pulse train 136 comprises the 
interference pulses 130, 132, and 134. 

The optical receiver 110 in one example comprises one or more photodiodes 138. In 
a fiirther example, the optical receiver 110 comprises a transimpedance amplifier 140. The 



6 NGC-144/000169-199 

optical receiver 110 in one example comprises a polarization diversity receiver system as 
disclosed in U.S. Patent No. 5,852,507 to Hall, issued May 10, 1996, and assigned to Litton 
Systems Inc. The optical receiver 110 receives the optical pulse train 136. The optical 
receiver 110 then creates one or more respective analog electrical signals that represent the 
5 interference pulses 130, 132, and 134 from the optical pulse train 136. For example, the 
optical receiver 110 converts a magnitude of power of the optical pulse train 136 to a voltage 
signal. 

The processor component 1 12 in one example comprises a digital signal processor. In 
a further example, the processor component 112 comprises an analog-to-digital converter 

10 component 142. The processor component 112 in one example comprises an instance of a 
computer-readable signal-bearing media 144, as described herein. The analog-to-digital 
converter component 142 converts the analog electrical signal from the optical receiver 110 
into a digital signal. The processor component 1 12 in one example serves to sense a change 
in the parameters by employing the time-varying amplitude variation of the interference 

15 pulses 130, 132, and 134 to calculate the time-varying phase angle cp. 

An illustrative description of exemplary operation of the apparatus 100 is presented, 
for explanatory purposes. The laser 102, the optical switch 104, and the phase modulator 106 
cooperate to create the one or more optical pulses 116. The sensor array 108 splits the optical 
pulse 116 into the optical pulses 118, 120, and 122. The sensors 124, 126, and 128 employ 

20 the parameters and the optical pulses 118, 120, and 122 to create the interference pulses 130, 
132, and 134. The sensor array 108 sends the interference pulses 130, 132, and 134 as the 
optical pulse train 136 to the optical receiver 1 10. 

The optical receiver 110 creates an analog electrical signal that represent the one or 
more interference pulses 130, 132, and 134. For example, the analog electrical signal is 

25 defined as s(t, M, P, <p): 
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where A is an average signal level, B is an interference term signal level, M is the 
modulation depth, Tpgc is the period of the phase generated carrier, p is the demodulation 
phase offset, and 9 is the phase angle. The phase angle of s(t, M, P, 9) comprises a first 



5 portion due to the phase generated carrier, M • sin 
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, and a second portion due to 



the parameter, 9, as will be understood by those skilled in the art. 

The analog-to-digital converter component 142 in one example converts the analog 
electrical signal from the optical receiver 110 into a digital signal that represents the 
interference pulse 130. The processor component 112 obtains only four samples of the 
10 interference pulse 130 from the digital signal. For example, the four samples comprise 
samples So, Si, S2, and S3. The processor component 112 obtains the four samples at time 
intervals At over a period Tg. The period Tg in one example is substantially equal to the 
period Tpgc of the phase generated carrier 114. The time interval At is equal to Tpgc / 4, For 
example, the four samples comprise samples at to, to + At, t© + 2At, to + 3 At. 

The processor component 112 employs one or more of the four samples to calculate 
an in-phase term 1. For example, the processor component 112 calculates the in-phase term I: 

I = (So + S2)-(S,+S3 ). 
The processor component 112 calculates a peak value Ip of the in-phase term 1. For example, 
the processor component 112 calculates the peak value Ip: 

Ip (M, /?) = 2 • B • (cos(M • sin 0) - cos(m • sin(^ + /?))). 
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The processor component 112 employs one or more of the four samples to calculate a 
quadrature term Q. In one example, the processor component 1 12 calculates the quadrature 
term Q: 

Q = _(So-S2). 

5 In another example, the processor component 1 12 calculates the quadrature term Q: 
Q=-2.(So-S,). 

The processor component 112 calculates a peak value Qp of the quadrature term Q. In 
one example where Q = -( So - Si ), the processor component 1 12 calculates the peak value 
Qp: 

10 Qp (M, >ff) = 2 • B • sin(M • sin p) 

In another example where Q = — 2 • (Sq — S2) , the processor component 112 calculates the 
peak value Qp: 

Q p (M, yff ) = 4 • B • sin(M • sin p) . 

The phase angle cp is given by: 

15 ^ = R(M,/?)-tan^, 

where a ratio R(M, P) is equal to Qp / Ip, as will be appreciated by those skilled in the art. 

The ratio R(M, P) in one example describes a ratio of principal axes of an elliptical 
Lissajous pattern for an operating point. The operating point in one example comprises a 
value of the modulation depth M and a value of the demodulation phase offset p. The 
20 processor component 112 in one example calculates the values for the modulation depth M 
and the demodulation phase offset P such that the ratio R(M, P) is substantially equal to one, 
to promote an increase in accuracy of a calculation of the phase angle 9. In one example 
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where Q = -( So - S2 ), the processor component 112 employs a value of 2.75 radians for the 
modulation depth M to calculate the phase offset P of the operating point, as described herein. 
In another example where Q = — 2-(So -83), the processor component 112 employs a 

value of 2.49 radians for the modulation depth M to calculate the phase offset (3 of the 
5 operating point, as described herein. 

The processor component 112 in one example employs the in-phase term I and the 
quadrature term Q to calculate the phase angle (p. For example, near an operating point, the 
ratio R(M, p) is substantially equal to unity, and the processor component 1 12 calculates the 
phase angle cp: 

10 (p == arctangent( Q / 1 ). 

The processor component 112 in one example employs the change in the phase angle (p 
between multiple instances of the interference pulses 130, 132, and 134 to determine the 
change in the parameters employed by the sensors 124, 126, and 128. 

Turning to FIG. 2-4, contour plots 202, 302, and 402 comprise a representation of the 

15 ratio R(M, p) for a previous design, a first implementation of the apparatus of FIG. 1, and a 
second implementation of the apparatus of FIG. 1, respectively. The contour lines represent 
the ratio Q / I versus the modulation depth M and the demodulation phase offset p. Contour 
line 204, 304, and 404 represent values for the modulation depth M and the demodulation 
phase offset P where the value of R(M, P) equals unity (100%). A modulation depth M is 

20 calculated to promote a decrease in sensitivity to the modulation depth M, as will be 
appreciated by those skilled in the art. For example, where a slope of the contour line 204, 
304, and/or 404 approaches zero, the ratio R(M, P) is least sensitive to the modulation depth 
M. 
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An accuracy of the calculation of the phase angle 9 comprises a relationship with a 
maximum amplitude of the quadrature term Q and the in-phase term I. For example, at an 
operating point, the ratio R(M, P) is equal to unity, and the peak value Qp and the peak value 
Ip are substantially equal. 

5 Referring to FIG. 2, the contour plot 202 comprises a representation of the ratio R(M, 

P) for a previous design that employs six samples. The contour line 204 comprises a value of 
"100%" for the ratio R(M, P). The contour Une 204 approaches zero at a modulation depth M 
of approximately 3.25 radians. A corresponding value of the demodulation phase offset p is 
0.2914 radians. A common operating point for the previous design comprises a modulation 

10 depth M of 3.14 radians and a demodulation phase offset P of 0.2908. Four of the six 
samples are used to calculate Q, and four of the six samples are used to calculate 1. The 
maximum amplitude for the quadrature term Q and the in-phase term I is equal to 4xB, where 
B is the interference term signal level. At the operating point, the amplitude of the quadrature 
term Q and the in-phase term is equal to: 

15 Qp(7i, 0.2908) = Ip(7i, 0.2908) = 3.24xB. 

Referring to FIG. 3, the contour plot 302 comprises an exemplary representation of 
the ratio R(M, P) for an exemplary implementation of the apparatus of FIG. 1 where Q = -( 
So - S2 ). The contour line 304 comprises a value of "100%" for the ratio R(M, P). The 
processor component 112 in one example calculates the modulation depth M approximately 

20 equal to 2.75 radians such that the slope of the contour line 304 approaches zero. The 
processor component 112 in a further example calculates a corresponding value of the 
demodulation phase offset p approximately equal to 0.5073 radians. At the operating point, 
Ip is substantially equal to Qp. For example, 

2 • B • (cos(M • sin yff)- cos(m • sin(^ + /?)))= 2 • B • sin(M • sin fi) . 
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Therefore, the maximum amplitude is equal to 2xB, as will be understood by those skilled in 
the art. At the operating point, Qp(2.75, 0.5073) = Ip(2.75, 0.5073) = l,97xB. 

Referring to FIG. 4, the contour plot 402 comprises an exemplary representation of 
the ratio R(M, p) for another exemplary implementation of the apparatus of FIG. 1 where Q = 
5 - 2 • (So - S2 ) . The contour line 404 comprises a value of "100%" for the ratio R(M, p). 

The processor component 112 in one example calculates the modulation depth M 
approximately equal to 2.49 radians such that the slope of the contour line 404 approaches 
zero. The processor component 1 12 in a further example calculates a corresponding value of 
the demodulation phase offset P approximately equal to 0.3218 radians. At the operating 
10 point, Ip is substantially equal to Qp. For example, 

2 - B - (cos(M • sin /?)- cos(m • sin(^ + >&)))= 4 • B • sin(M • sin >5) . 

Therefore, the maximum amplitude is equal to 4xB, as will be understood by those skilled in 
the art. At the operating point, Qp(2.49, 0.3218) = Ip(2.49, 0.3218) - 2.83xB. 

Referring to FIGS. 2-4, the modulation depth M of the respective operating point 

15 comprises a relationship with a drive voltage of the phase generated carrier. A lower drive 
voltage promotes a decrease in power consumption of the phase generated carrier, as will be 
appreciated by those skilled in the art. The modulation depth of the implementation of FIG. 3 
is 88% of the modulation depth of the previous design of FIG. 2. The modulation depth of 
the implementation of FIG. 4 is 79% of the modulation depth of the previous design of FIG. 

20 2. 

Turning to FIG. 5-7, contour plots 502, 602, and 702 comprise exemplary 
representations of the contour plots 202, 302, and 402, respectively, with levels of total 
harmonic distortion. Contour lines 504, 604, and 704 comprise contour lines where R(M, p) 
is equal to "99%." Contour lines 506, 606, and 706 comprise contour lines where R(M, p) is 
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equal to "99.5%." Contour lines 508, 608, and 708 comprise contour lines where R(M, p) is 
equal to "100%." Contour lines 510, 610, and 710 comprise contour lines where R(M, p) is 
equal to "100.5%." Contour lines 512, 612, and 712 comprise contour lines where R(M, P) is 
equal to "101%." Total harmonic distortion is evaluated using a Hanning filter on the phase 
5 angle cp in the calculation of (p = arctangent (Q / I). A tone with a peak-to-peak amplitude 
equal to n is used to evaluate the total harmonic distortion. For example: 



where ftone comprises a frequency of the tone, a comprises a tone offset phase shift, and (pdc 
comprises a non-alternating phase angle of the tone. The total harmonic distortion of the 2"** 

10 through 10^*^ harmonics is evaluated at the operating point for the calculation of cp of FIGS. 2- 
4 where the total harmonic distortion is at a minimum. 

Referring to FIG. 5, contour plot 502 comprises a representation of the ratio R(M, P) 
of FIG, 2. The frequency fpgc is equal to 24.576 kHz and ftone is equal to 96 Hz. 4096 
samples were taken over a period of 166.67 milliseconds. Values for total harmonic 

15 distortion for the contour lines 504, 506, 508, 510, and 512 comprise -49 dB, -54 dB, -59 dB, 
-54 dB, and -49 dB, respectively. 

Referring to FIG. 6, contour plot 602 comprises a representation of the ratio R(M, P) 
of FIG. 3. The frequency fpgc is equal to 36.864 kHz and ftone is equal to 99 Hz. 4096 
samples were taken over a period of 111.11 milliseconds. Values for total harmonic 

20 distortion for the contour lines 604, 606, 608, 610, and 612 comprise -49 dB, -53 dB, -57 dB, 
-53 dB, and -49 dB, respectively. 

Referring to FIG. 7, contour plot 702 comprises a representation of the ratio R(M, P) 
of FIG. 4, The frequency fpgc is equal to 36.864 kHz and ftone is equal to 99 Hz, 4096 
samples were taken over a period of 111.11 milliseconds. Values for total harmonic 
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distortion for the contour lines 704, 706, 708, 710, and 712 comprise -49 dB, -53 dB, -56 dB, 
-53 dB, and -49 dB, respectively. 

Referring to FIGS 5-7, the accuracy of the calculation of the phase angle <p further 
comprises a relationship with a minimum horizontal shift in the modulation depth M and a 
5 minimum vertical shift in the demodulation phase offset P between the contour 100 and the 
contour 100.5. The calculation of FIG. 5 comprises a shift in modulation depth M of 0,095 
radians and a shift in demodulation phase offset P of 0.0015 radians. The calculation of FIG. 
6 comprises a shift in modulation depth M of 0.16 radians and a shift in demodulation phase 
offset P of 0.0013 radians. The calculation of FIG. 7 comprises a shift in modulation depth 

10 M of 0. 1 6 radians and a shift in the demodulation phase offset P of 0.00 1 3 radians. 

The accuracy of the calculation of the phase angle cp for the previous design of FIGS. 
2 and 5 and the exemplary implementations of FIGS. 3, 4, 6, and 7 comprise a similar 
accuracy. The processor component 112 in one example employs only four samples as 
opposed to six samples to promote an increase in speed of calculation of the phase angle 9. 

15 For example, the processor component 112 calculates the phase angle 9 50% faster than the 
previous design. The increase in sampling rate allows for an increase in frequency fpgc of the 
phase generated carrier, as will be appreciated by those skilled in the art. 

The apparatus 100 in one example employs one or more computer-readable 
signal-bearing media. An exemplary computer-readable signal-bearing media for the 

20 apparatus 100 comprises the recordable data storage media 144 of the processor component 
112. For example, the computer-readable signal-bearing media for the apparatus 100 
comprises one or more of a magnetic, electrical, optical, biological, and atomic data storage 
media. In one example, the computer-readable signal-bearing media comprises a modulated 
carrier signal transmitted over a network comprising or coupled with the apparatus 100, for 
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instance, one or more of a telephone network, a local area network ("LAN*'), the internet, and 
a wireless network. 

The steps or operations described herein are just exemplary. There may be many 
variations to these steps or operations without departing from the spirit of the invention. For 
5 instance, the steps may be performed in a differing order, or steps may be added, deleted, or 
modified. 

Although exemplary implementations of the invention have been depicted and 
described in detail herein, it will be apparent to those skilled in the relevant art that various 
modifications, additions, substitutions, and the like can be made without departing from the 
10 spirit of the invention and these are therefore considered to be within the scope of the 
invention as defined in the following claims. 



